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Abstract 

We  present  a  first-principles  study  to  understand  the  phenomena  of  interlayer  exchange  cou¬ 
pling  in  Fe/Nb  multilayers  using  the  linearized- muffin-tin-orbitals  method  within  the  generalized 
gradient  approximation.  We  find  that  the  exchange  coupling  oscillates  with  both  short  and  long 
periodicities,  which  have  been  examined  in  terms  of  the  Ruderman-Kittel-Kasuya-Yosida  (RKKY) 
model  as  well  as  the  quantum  well  (QW)  model.  We  have  investigated  the  behavior  of  the  ex¬ 
change  coupling  by  artificially  varying  moments  of  Fe  atoms  in  ferromagnetic  layers.  For  small 
moment  of  Fe,  the  coupling  shows  bilinearity  in  the  magnetic  moments  implying  its  RKKY  charac¬ 
ter.  However,  at  higher  moments  close  to  the  bulk  Fe,  the  saturation  of  long-period  oscillations  is 
in  accordance  with  the  QW  model.  Quantum-well  dispersions  around  the  Fermi  level  demonstrate 
that  the  majority-spin  bands  contribute  largely  to  the  formation  of  quantum-well  states,  which 
we  analyze  quantitatively  by  making  use  of  the  phase  accumulation  model.  Our  analysis  indicates 
that  the  quantum  well  model  gives  a  better  description  of  the  oscillatory  behavior  of  the  exchange 
coupling  in  Fe/Nb  multilayers. 
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I.  INTRODUCTION 


Damped  long-range  oscillation  of  the  interlayer  exchange  coupling  (lEC)  as  a  function 
of  the  spacer  thickness  is  a  well  known  phenomenon  in  magnetic  multilayers [1-3].  Several 
approaches  have  been  proposed  over  the  years  in  order  to  explain  the  oscillatory  behavior 
of  the  lEC.  Among  these,  the  two  prominent  models  are  (i)  the  Ruderman-Kittel-Kasuya- 
Yosida  (RKKY)  model[4]  and  (ii)  the  quantum  well  (QW)  model[5].  The  RKKY  interaction 
stems  from  the  spin-polarization  of  the  intervening  conduction  electrons  in  the  spacer  layer. 
Bruno  and  Chappert[6]  have  shown  that  the  exchange  coupling  within  the  RKKY  theory  is 
related  to  the  topological  properties  of  the  Fermi  surface  of  the  spacer  material.  Since  the 
magnetic  atoms  in  multilayers  are  more  immersed  in  the  spin  sea  of  the  conduction  electrons 
of  the  spacer  layer  near  the  interface,  the  RKKY  interaction  has  the  strongest  effect  at  the 
interfaces [7].  However,  the  RKKY  approach  is  not  much  effective[4,  6]  in  describing  the 
correct  amplitude  and  phase  of  the  coupling  oscillations  essentially  because  the  interaction 
between  the  ferromagnetic  layers  and  the  conduction  electrons  is  not  well  captured  in  the 
RKKY  approach.  In  the  Quantum  Well  model,  the  coupling  arises  due  to  the  spin-dependent 
confinement  [8]  of  electrons  inside  the  spacer  medium  as  the  size  of  the  multilayer  system 
is  reduced  to  the  nanometer  range.  In  the  quantum  well  picture,  each  layer  thickness  of 
magnetic  as  well  as  nonmagnetic  (NM)  kind  in  the  entire  multilayer  stack  contributes  sig¬ 
nificantly  to  the  coupling  strength  implying  that  the  lEC  is  not  a  sheer  interfacial  effect  [9]. 
The  QW  states  are  formed  near  the  Fermi  level  when  spin-polarized  bands  are  shifted  away 
from  the  Fermi  surface  because  of  the  strong  magnetization  on  both  sides  of  the  spacer 
medium.  Such  states,  which  can  be  observed  directly  by  photoemission  measurements,  shift 
in  energy  with  the  spacer  thickness  and  become  closely  spaced  when  the  spacer  layer  ap¬ 
pears  to  be  sufficiently  thick.  Though  both  the  models  arrive  at  the  same  period  of  coupling 
oscillations  due  to  their  origin  in  the  shape  of  the  spacer  Fermi  surface,  they  differ  in  deter¬ 
mining  the  coupling  strength [10].  It  is  because  the  RKKY  model  originates  from  the  second 
order  perturbation  theory  whereas  the  QW  model  does  not  make  any  assumption  about  the 
strength  of  the  interlayer  interaction.  Further,  QW  theory  predicts  additional  possibilities 
like  oscillating  density  of  states,  quantum  well  dispersions  etc.,  especially  when  the  spacer 
layer  forms  a  multisheet  Fermi  surface  (FS).  Such  is  the  case  in  Fe/Nb  multilayers,  where 
Nb  spacer  has  three  sheets  of  Fermi  surface  in  the  (100)  plane[ll]. 
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Although  many  computational  works  exist  already  in  the  literature[2,  4,  12]  on  the  ex¬ 
change  coupling  phenomena  in  systems  such  as  Fe/Cr,  Fe/Au,  Co/Cu  etc,  where  the  spacer 
layers  are  in  general  NM  transition  metals,  not  much  theoretical  understanding  has  been 
gained  so  far  for  Fe/Nb  multilayer  system.  The  first  experimental  study  on  the  sputtered 
Fe/Nb  superlattices  as  carried  out  by  Mattson  et  a/ [13]  reported  a  weak  coupling  with  an 
oscillation  period  of  about  9  A  at  room  temperature.  However,  some  ab  initio  band  struc¬ 
ture  studies[14,  15]  on  Fe/Nb  multilayers  demonstrated  an  oscillation  period  of  4. 6-6.0  A. 
Neutron  reflectometry  data  of  Rehm  et  a/ [16]  for  Fe/Nb  multilayers  have  earlier  suggested 
an  oscillatory  RKKY  kind  of  coupling  for  small  Nb  layers.  In  another  experimental  de¬ 
velopment,  Klose  et  al[17]  have  shown  that  hydrogen  charging  can  modify  the  magnetic 
coupling  in  these  heterostructures  through  the  alteration  of  the  electronic  structure  of  the 
Nb  interlayer. 

Fe/Nb  multilayers  seem  to  be  interesting  systems  to  study  the  interlayer  exchange  cou¬ 
pling  in  the  sense  that  Fe  is  a  strong  transition-metal  ferromagnet  (FM)  while  the  spacer 
layer  can  be  a  superconductor  (SC)  at  low  temperatures.  Thus  apart  from  the  phenomena 
of  oscillatory  exchange  coupling  between  FM  layers,  Fe/Nb  superlattices  have  another  inter¬ 
esting  phenomenon  known  as  the  ’’proximity  effect”,  which  paves  the  way  for  new  sources 
of  magneto-resistance  with  potential  applications  in  magneto-electronics[18].  The  presence 
of  the  internal  magnetic  field  in  Fe  layers  weakens  the  phenomena  of  superconductivity  due 
to  the  breaking  of  Cooper  pairs.  Among  various  FM/SC  hetero-structures,  where  supercon¬ 
ductivity  gets  induced  in  the  ferromagnet  by  bringing  it  in  contact  with  the  superconductor, 
Fe/Nb  multilayers  have  been  studied  experimentally[16,  19,  20]  to  an  appreciable  extent. 
One  such  observation  by  Miihge  et  a/ [20]  suggests  that  the  thickness  of  the  SC  interface 
mainly  determines  the  actual  shape  of  versus  dpe  (thickness  of  the  FM  layer)  curve  in 
Fe/Nb  superlattices.  Since  Nb  is  the  highest  known  element,  its  interplay  with  Fe  con¬ 
tinues  to  remain  an  active  field  in  the  study  of  such  heterostructures. 

In  the  present  work,  we  elucidate  the  coupling  phenomena  in  Fe/Nb  multilayers  within  the 
domain  of  the  density  functional  theory[21].  Previous  calculation  of  Shukla  and  Prasad[15] 
reported  the  lEC  for  the  Fe/Nb  multilayer  system  up  to  7  monolayers  of  Nb  spacer  sand¬ 
wiched  between  two  Fe  layers.  However,  it  could  not  shed  light  on  the  detailed  analysis  of 
the  coupling  phenomena  due  mainly  to  the  lack  of  relatively  large  spacer  thickness,  which  is 
necessary  for  such  kind  of  study.  This  work  thus  involves  FesNb^  (m=1..16)  system  in  order 
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to  understand  the  coupling  behavior  in  Fe/Nb(001)  multilayers.  Our  calculation  shows  that 
QW  states  are  indeed  formed  in  such  heterostructures  and  the  QW  model  gives  a  better 
description  of  the  oscillatory  exchange  coupling  in  Fe/Nb  multilayers. 

The  organization  of  the  paper  is  as  follows.  In  Sec.  II,  we  briefly  outline  the  computational 
procedure  adopted  in  the  present  study.  Sec.  Ill,  under  several  subsections,  deals  mainly 
with  our  results  that  include  simultaneous  discussions.  We  finally  sum  up  our  observations 
in  Sec.  IV. 

II.  METHODOLOGY 

Total  energy  calculations  have  been  carried  out  for  FesNb^  (m=1..16)  multilayers  in  the 
framework  of  the  linearized-muffin-tin-orbitals  (LMTO)  method[22-24]  within  the  tight- 
binding  representation.  The  atomic-sphere  approximation  (ASA)  is  used  for  the  potentials 
determined  self  consistently  in  the  generalized  gradient  approximation  (GGA)[25]  of  the 
density  functional  theory [21].  Various  tetragonal  supercells  are  constructed  out  of  Fe  and 
Nb  monolayers  (ML)  where  the  bcc  Fe  layers  are  stacked  along  the  [001]  growth  direction 
(see  Ref.  15)  in  ferromagnetic  (FM)  as  well  as  antiferromagnetic  (AFM)  orientations.  The 
interlayer  exchange  coupling  (lEC  ),  denoted  by  J(m),  corresponds  to  the  energy  difference 
between  FM  and  AFM  configurations  per  unit  cell  structure  so  that 

J (m)  =  Eli{m)  -  Et'^t(m),  (1) 

where  m  is  the  number  of  spacer  layers.  We  calculate  the  total  energies  of  all  occupied  states 
and  minimize  it  between  the  FM  and  AFM  configurations  in  the  self-consistent  fashion  for 
each  Nb  thickness.  The  average  lattice  parameter  with  reduced  lattice  mismatch  is  taken [15] 
as  3.067  A  for  the  present  heterostructures.  Linear  tetrahedron  method  has  been  used  for  the 
Brillouin  zone  (BZ)  integration  with  a  maximum  of  840  /c-points  in  the  irreducible  wedge  of 
the  surface  BZ.  We  use  the  same  unit  cell  for  the  FM  and  AFM  structures  to  obtain  reliable 
energy  differences  between  ferromagnetic  and  antiferromagnetic  ordering. 
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III.  RESULTS  AND  DISCUSSION 


A.  Interlayer  coupling  oscillations 

Initially,  we  compute  the  lEC  as  a  function  of  the  Nb  spacer  thickness  using  Eq.  (1) 
for  FeNb^(m  =  1..7)  system.  The  results  for  J{m)  are  in  good  agreement  with  our  earlier 
calculation [15]  using  the  FP-LAPW  method.  This  gives  us  confidence  in  carrying  out  the 
present  LMTO-ASA  based  calculations,  which  provide  a  reasonable  estimate  of  magnetic 
moments  and  energy  differences [26].  The  lEC  for  Fe3Nbm(m  =  1..16)  configuration  is  then 
computed  using  Eq.  (1),  which  is  shown  in  Fig.  1.  Rapid  oscillations  are  observed  up  to 
9  monolayers  (ML)  of  Nb  thickness,  after  which  the  oscillatory  exchange  coupling  becomes 
appreciably  weak.  The  coupling  changes  from  ferromagnetic  to  antiferromagnetic  configu¬ 
ration  at  about  2,  7  and  10  monolayers  of  Nb  spacer,  resulting  in  oscillation  periods  of  7.7 
and  4.6  A  respectively.  Note  that  1  ML  corresponds  to  the  interplanar  thickness  of  1.5335 

A. 

To  better  understand  the  oscillation  periods  and  the  coupling  phenomena  in  Fe/Nb  mul¬ 
tilayers,  we  fit  the  calculated  variation  of  J{m)  with  m  to  the  following  asymptotic  form[8] 

2  4 

J(m)  =  E  Aksin  {qkfn  -F  (f)k)  Im^  -F  ^  Aisin  {qim  -F  /m^  (2) 

k=l  1=3 

where  Aa;(/)’s  describe  the  amplitudes;  qk(i)^s  yield  the  periodicities  Tk(i){=2'K / qk[i)ys]  and 
(t>k(iys  the  phases  of  the  A:(/)-th  mode  of  oscillations.  Our  ab  initio  data  for  J(m)  in  Eq. 
(1)  yields  a  reasonably  good  fit  to  Eq.  (2)  with  the  following  four  periods:  Ti  =  4.14  ML 
6.3  A),  Ta  =  5.05  ML  (~  7.7  A),  Tg  =  2.86  ML  (~  4.4  A),  T4  =  20.28  ML  (~  31.1  A). 
These  values  fall  well  within  the  previous  results[13,  15]  for  the  FeNb^  heterostructures. 
Multiple  periodicities  arise  due  to  existence  of  the  multisheet  Fermi  surface  in  the  Nb  spacer 
layer.  The  Nb  thickness  periodicity  of  7.7  A  comes  close  to  the  experimental  value  of  9.0 
A.  This  discrepancy  is  perhaps  due  to  pre-asymptotic  effects  and  the  difficulty  of  including 
the  true  lattice  structure  in  our  calculations.  Since  in  the  present  study,  the  total-energy 
calculations  of  the  lEC  are  limited  by  16  ML  of  spacer  thickness,  which  might  still  be  in  the 
pre-asymptotic  region,  a  slight  discrepancy  between  the  experimental  result  and  total-energy 
calculation  is  expected.  However,  the  4.4  and  6.3  A  periods  are  in  fairly  good  agreement 
with  the  available  ab  initio  data[13-15,  27]. The  well-fitted  curve  of  Fig.  1  suggests  that 
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the  interlayer  exchange  coupling  in  Fe/Nb  multilayers  has  significant  1/w?  dependence  in 
addition  to  the  conventional  1/m^  dependence  given  by  the  RKKY  theory. 

If  we  analyze  the  four  periodicities  that  are  obtained  upon  a  reasonably  good  fitting  of 
Eq.  (2),  the  existence  of  higher  harmonics  is  observed  in  the  coupling  function.  These 
harmonics  (denoted  by  n)  add  terms  of  the  form  sin{2'Knm/T  -\- phase)  with  i  =  0, 

1...,  n  >  2,  and  T  being  the  fundamental  period.  This  way  the  fourth  harmonic  (n  =  5) 
results  in  m~'^ sin{56m  +  phase)  with  6  =  27^  jT.  Now,  if  we  express  the  4.14  ML  period  Ti 
as  1/Ti  =  1/5  +  5/27r,  the  effective  period[8]  is  obtained  as  T^ff  =  27r/55  =  Ti/(5  —  Ti)  = 
4.81  ML.  This  value  is,  however,  a  bit  lower  than  the  fit  value  of  5.05  ML,  which  corresponds 
to  T2.  The  slight  discrepancy  may  be  attributed  to  the  uncertainties  in  the  fit  by  Eq.  (2). 
Expressing  similarly  the  2.86  ML  period  T3  as  I/T3  =  1/3  +  5/27r,  the  second  harmonic 
(n  =  3)  becomes  m~^sin{35m  +  phase)  ^  which  yields  an  effective  period  =  2'k/35  = 
^3/(3  —  Ts)  =  20.43  ML.  This  value  is  also  quite  close  to  the  fit  value  of  T4  (=  20.28  ML) 
with  a  discrepancy  of  about  0.7%.  Both  the  long  periods  T2  and  T4  thus  turn  out  to  be  the 
’’Vernier”  periods  of  the  respective  short  period  oscillations  of  Ti  and  T3.  As  pointed  out 
by  Schilfgaarde  and  Harrison[8],  higher  harmonics  have  significant  presence  in  the  quantum 
well  limit  where  the  RKKY  description  does  not  hold  good. 

B.  Fermi  surface  and  the  RKKY  periods 

In  the  RKKY  approach,  the  oscillatory  periods  of  the  interlayer  exchange  coupling  are 
uniquely  determined  by  the  stationary  spanning  vectors  of  the  bulk  Fermi  surface  of  the 
spacer  material.  Several  spanning  vectors  in  the  FS  give  rise  to  multiple  periodicities,  as  we 
have  already  come  across  in  the  previous  section. 

In  Fe/Nb  multilayers,  the  spacer  layer  Nb  has  five  conduction  electrons  per  atom  that  fill 
the  first  Brillouin  zone  completely  while  the  second  and  third  Brillouin  zones  partially [28]. 
The  second  zone  is  a  closed  ” octahedron”  (OCT),  which  contains  a  hole  sheet  centered  at  F. 
However,  the  third  zone  has  two  sheets.  One  sheet  contains  an  open  surface  of  holes,  referred 
to  as  ”  jungle  gym”  (JG),  which  extends  from  F  to  iL  points  along  the  [100]  direction.  The 
other  sheet  is  a  set  of  distorted  hole  ”  ellipsoids”  (ELL)  centered  at  N  points.  The  OCT 
and  JG  sheets  contact  at  three  points[29]  in  the  (100)  and  (110)  symmetry  planes  of  the 
Fermi  surface  of  Nb. 
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Simple  square  lattice  planes  of  Fe  and  Nb  with  primitive  translations  of  ^[100]  and 
^[010]  are  stacked  along  the  [001]  growth  direction  to  form  Fe/Nb  multilayers.  Since  the 
translational  symmetry  in  the  growth  direction  is  broken,  multilayers  behave  like  quasi-two- 
dimensional  systems,  which  are  periodic  only  in  two  dimensions[30].  Hence,  the  Brillouin 
zone  needs  to  be  constructed  in  2D  to  deal  with  the  in-plane  coordinates[31].  Fig.  2  shows 
the  cross  section  of  the  spacer  layer  Fermi  surface  of  Nb  on  the  basis  of  our  self-consistent 
calculations  in  the  central  (100)  plane.  The  line  Qi,  Q2  and  as  shown  in  Fig.  2,  represent 
the  spanning  vectors  for  the  [100]  crystalline  orientation;  Qi  spans  the  F-centered  octahedron 
along  the  [100]  direction,  while  Q2  and  span  the  outer  and  the  inner  ellipses  respectively 
along  the  [100]  direction.  We  find  that  the  RKKY  periods  of  4.1,  6.4  and  7.5  A,  as  predicted 
from  the  Fermi  surface  spanning  vectors,  turn  out  to  be  in  good  agreement  with  the  interlayer 
coupling  periods  of  4.4,  6.3  and  7.7  A,  as  obtained  by  fitting  the  self-consistent  results  to 
Eq.  (2).  However,  the  Fermi  surface  topology  of  Nb  is  not  able  to  predict  the  periodicity 
of  31.1  A.  As  we  know  from  the  preceding  section,  the  long-wavelength  coupling  period  of 
31.1  A  originates  from  the  higher  order  terms  in  the  coupling  function  of  Eq.  (2),  which  is 
also  a  higher  harmonics  of  the  short-wavelength  period  of  4.4  A  . 

C.  Magnetization 

To  examine  how  the  Fe  magnetic  moment  at  the  interface  behaves  as  a  function  of  the 
intervening  layer  thickness  in  Fe/Nb  multilayers,  we  plot  the  magnetic  moment  of  the  in¬ 
terface  monolayer  of  Fe  vs.  the  spacer  layer  thickness  of  Nb  in  Fig.  3.  There  is  a  good 
agreement  between  our  calculated  results  and  the  experimental  observations  as  shown  in 
the  inset.  We  find  that  the  Fe  magnetic  moment  gets  saturated  in  the  asymptotic  region 
with  reduction  of  about  25  %  of  the  bulk  value.  This  reduction  in  the  measured  data  is 
about  40  %.  Following  the  works  of  Holmstrom  et  a/ [32],  this  kind  of  discrepancy  may  be 
attributed  primarily  to  the  interface  alloying. 

The  induced  magnetic  moments  in  the  Nb  spacer  layer  for  a  ferromagnetically  ordered 
FesNbie  superlattices  are  shown  in  Fig.  4(a).  On  the  other  hand.  Fig.  4(b)  displays 
the  induced  polarization  for  the  antiferromagnetic  configuration.  Fe  layers  are  at  positions 
marked  by  0  and  17  (not  shown  in  the  figure).  The  induced  magnetic  moment  (in  absolute 
values)  in  Nb  spacer  layer  decreases  from  about  0.28  [Ib  at  the  interface  to  about  0.02  [Ib 
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further  apart  from  the  interface.  For  an  Nb  layer  of  16  atomic  planes,  the  calculated  induced 
moment  at  the  center  of  the  Nb  layer  is  about  10“^  The  period  of  oscillation  of  the 
induced  moment  in  both  the  magnetic  configurations  turns  out  to  be  about  4.6  A,  which  is 
in  fairly  good  agreement  with  the  coupling  period  of  4.4  A,  as  determined  by  Eq.  (2).  The 
bias  in  the  magnetic  ordering  of  the  induced  moment  (see  Fig.  4),  which  occurs  only  in  the 
ferromagnetic  configuration,  may  be  regarded  as  being  due  to  the  onset  of  non-RKKY  terms 
in  the  coupling  function  when  the  Nb  layer  gets  appreciably  thicker.  Mathon  et  a/ [33]  have 
already  shown  analytically  the  presence  of  such  non-RKKY  terms  for  Co/Cu(001)  system 
using  the  stationary  phase  approximation. 

However,  to  get  a  better  understanding  of  the  coupling  behavior,  one  needs  to  examine 
the  influence  of  Fe  magnetization  on  the  coupling  strength.  For  this,  we  adopt  a  procedure 
similar  to  that  of  Schilfgaarde  and  Harrison[8].  A  trial  density  is  constructed  out  of  the 
charge  densities  of  self-consistently  calculated  bulk  ferromagnetic  Fe  and  paramagnetic  Nb  in 
their  respective  atom-centered  spheres.  The  charge  density  of  Fe  is  constructed  as  follows[8] 

n,.M=«»M±adk)^,  (3) 

where  denotes  the  density  of  bulk  paramagnetic  Fe  while  and  nf  represent  the  spin 
densities  of  the  majority  and  minority  spins  in  ferromagnetic  Fe  and  cr  is  a  parameter  ranging 
from  0  to  1.  The  magnetic  moment  of  Fe  atom  will  be  proportional  to  a,  taking  full  value  at 
q;  =  1  and  0  value  when  a  =  0.  Since  we  are  interested  in  obtaining  J  as  a  function  of  the 
Fe  moment,  we  have  used  the  trial  densities  in  our  frozen-potential  calculations  pertaining 
to  relevant  supercells.  To  have  a  preliminary  idea  of  how  the  coupling  and  the  moments 
are  inter-related  in  Fe/Nb  systems,  we  construct  a  16-atom  Nb  supercell  with  two  Fe  atoms 
substituted  such  that  one  Fe  sits  at  (0  0  0)  while  the  other  at  3a(l  1  l)/2.  For  parallel  and 
antiparallel  alignment  of  two  inequivalent  Fe  atoms,  the  energy  difference  is  calculated [34] 
as  a  function  of  the  Fe  moment  parametrized  by  a.  According  to  the  RKKY  theory,  this 
energy  difference  should  vary  as  (see  Fig.  5). 

Several  superlattices  out  of  FeaNb^  (m  =  1-16)  multilayer  configuration  are  con¬ 
structed  subsequently  to  calculate  the  energy  difference  Ja{fn)  =  -^'[FegNbmFegNb^]  - 
EpegNb^FegNbrn],  as  function  of  a.  We  determine  the  amplitudes  A;t(/)’s  for  each  ol  by 
fitting  JaiiTi)  to  the  functional  form  as  given  by  Eq.  (2).  This  way,  we  find  two  short-period 
amplitudes  and  two  long-period  amplitudes,  which  are  illustrated  in  Fig  5.  In  Fe/Cr 
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multilayers,  Schilfgaarde  and  Harrison [8]  found  the  coupling  amplitudes  (for  small  a)  of 
both  short  and  long  periods  to  increase  as  while  at  larger  a,  short  period  amplitude 
continued  its  rise,  though  the  long  period  saturated  completely.  Here  we  find  that  the 
initial  dependence  of  the  exchange  coupling,  as  assumed  in  the  RKKY  theory,  at  small 
OL  values  is  followed  by  two  distinct  features  at  higher  ol  values  viz. 

(i)  stronger  dependence  on  the  Fe  moment  as  lEC  shoots  up  for  short  period  oscillations 
and 

(ii)  complete  independence  of  the  Fe  moment  as  lEC  gets  saturated  for  long  period 
oscillations. 

As  a  function  of  magnetization  {/j,),  the  interlayer  exchange  energy  can  be  described  as[35] 

jjj,  =  JlitJ'lfJ-2)  +  (4) 

where  Ji  and  J2  are  respectively  the  bilinear  and  biquadratic  coupling  constants,  while  /ii 
and  )U2  denote  the  magnetizations  of  the  adjacent  ferromagnetic  layers.  In  order  to  fit  the 
calculated  results  for  the  16-atom  supercell  data,  we  need  to  expand  the  above  expression 
in  the  following  form 

i/q,  =  J\0?  -|-  J20i^  “1“  (5) 

where  (i  has  been  parametrized  by  a  and  J3  is  an  additional  triquadratic  term  appearing 
in  the  exchange  coupling.  Initially,  we  calculate  the  coupling  energy,  as  a  function  of  a,  for 
the  16-atom  Nb  supercell  with  only  two  Fe  atoms  substituted.  These  are  represented  by 
asterisks  in  Fig.  5.  While  fitting  the  data  we  find  that  the  coupling  energy  has  biquadratic 
and  triquadratic  terms  in  addition  to  bilinear  terms  in  a.  In  doing  so,  the  coupling  constants 
turn  out  to  be  in  the  ratio  of  about  1:2:4.  It  may  be  noted  that  the  RKKY  theory  assumes 
only  the  bilinear  terms  in  a.  The  linear  part  of  the  fitted  curve  in  Fig.  5  thus  shows  the 
RKKY  kind  of  coupling  in  the  region  of  small  magnetic  moment  of  Fe.  We  then  calculate 
the  interlayer  exchange  coupling  for  the  FesNb^  (m  =  1-16)  multilayer  configuration  as  a 
function  of  a.  Short  period  amplitudes  of  the  FesNb^  (m  =  1-16)  multilayer  configuration 
are  found  to  follow  the  above  phenomenological  expression  for  the  interlayer  interaction, 
but  the  long  period  amplitudes  show  a  saturation  behavior  at  large  a.  For  a  complete 
confinement,  the  coupling  becomes  independent  of  the  size  of  the  magnetic  moment [36],  as 
demonstrated  by  the  long  period  amplitudes.  The  saturation  of  the  long-period  amplitudes 
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with  the  ferromagnetic  Fe  moment  thus  favors  a  QW  description  of  the  exchange  coupling. 
Because  of  the  incomplete  confinement  of  states  in  Nb,  Fe/Nb  multilayers  exhibit  a  partial 
transition  to  the  QW  character  from  the  RKKY  description,  which  is  valid  at  small  moments 
[Fig.  5]. 


D.  Energy  bands,  density  of  states  and  the  quantum  well  model 

To  study  the  nature  of  the  electronic  states  in  Fe/Nb  multilayers,  we  resort  to  the  bulk 
band  structures  (see  Fig.  6)  of  the  ferromagnetic  Fe  (for  both  spins)  and  the  paramagnetic 
Nb  along  the  [100]  as  well  as  the  [110]  direction.  These  directions  are  important [11]  for 
multiple  scattering,  since  the  respective  A2,  Si  and  S2  bands  cross  the  Fermi  level  (E^’)  in 
the  spacer  medium  contributing  to  the  oscillations  in  the  interlayer  coupling. 

The  Fe  minority  F12-A2-H12  band  closely  resembles  the  corresponding  Nb  band  indicating 
that  the  Fe/Nb  interface  will  be  more  transparent  to  the  minority-spin  electrons  than  to  the 
majority-spin  electrons.  The  vanishing  overlap  between  the  majority  spin  bands  Fi2-A2-i?i2 
in  Nb  and  Fe  creates  spin-dependent  gaps  at  the  Fe/Nb  interface  resulting  in  energy  barriers 
that  confine  electrons  to  the  intervening  Nb  layer.  The  electrons  in  Nb  band  Fi2-A2(t)--ffi2 
can  not  pass  into  the  Fe  majority  band  since  there  are  no  available  states  of  A2(t)  symmetry 
from  Fi2  (Fet)  to  F12  (Nb).  As  a  result,  spin-up  electrons  with  energies  from  (E^?  -  0.14) 
eV  to  (Ej^  -|-  3.15)  eV  experience  multiple  reflections  at  the  interface  and  get  confined  to 
the  Nb  spacer.  Similar  phenomena  happen  along  the  [110]  direction.  The  electrons  in  the 
Nb  band  F25'-Si(j")-W'  can  not  pass  into  the  Fe  majority  band  in  the  energy  range  from 
(Eji’  -  0.94)  eV  to  (E^?  -|-  2.38)  eV.  Similarly,  the  Nb  electrons  in  F25'-S2(t)-A^2  band  can  not 
pass  into  the  Fe  majority  band  in  the  energy  range  from  (E^’  -  2.29)  eV  to  (Ej?  -|-  0.56)  eV. 
The  confinement  of  electrons  thus  occurs  in  the  spacer  layer  resulting  in  the  formation  of 
quantum  well  states.  As  Fig.  6  shows,  bcc  Fe  has  mostly  minority-spin  states  at  the  Fermi 
level  and  this  gives  rise  to  majority-spin  quantum  well  states  in  Nb. 

Fig.  7  displays  how  the  density  of  states  at  Ep  oscillates  with  the  spacer  layer  thickness, 
which  is  a  characteristic  feature  of  the  quantum  well  model.  According  to  the  QW  model, 
the  oscillatory  behavior  of  the  magnetic  coupling  in  Fe/Nb  multilayers  as  demonstrated  by 
Fig.  1  can  be  traced  back  to  the  oscillations  of  the  density  of  states  (DOS)  at  the  Fermi 
level[5,  27].  During  the  formation  of  QW  states,  Fe/Nb  interfaces  act  as  electron  mirrors 
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that  induce  standing  waves  in  the  Nb  spacer  medium.  As  Fig.  7  suggests,  the  first  two 
maxima  correspond  to  the  oscillation  period  of  about  4.6  A.  The  periodicity  then  enhances 
to  about  6.1  A.  We  notice  that  the  QW  periodicities  of  4.6  and  6.1  A  are  in  fairly  good 
agreement  with  the  interlayer  coupling  periods  of  4.4  and  6.3  A  as  obtained  from  the  self- 
consistent  results  for  FeaNb^  (m=1..16)  multilayers  (see  Sec.  Ill- A).  It  implies  that  there 
is  a  reasonable  connection  between  the  existence  of  the  spin-polarized  QW  states  and  the 
manifestation  of  magnetic  coupling  in  Fe/Nb  multilayers. 

Oscillations  in  the  density  of  states  at  the  Fermi  level  can  also  be  explained  on  the  basis 
of  the  bulk  band  structures  of  Fe  and  Nb,  assuming  the  QW  model.  We  know  from  the 
previous  discussion  that  the  majority-spin  energy  bands  of  Fe  along  the  [100]  and  the  [110] 
directions  provide  the  confinement  to  Nb  electrons  in  forming  the  QW  states  inside  the 
spacer  medium.  Fig.  6  shows  that  A2  and  Si  bands  lead  to  the  formation  of  QW  states 
around  Ep  inside  the  Nb  layer.  It  is  because  there  are  no  propagating  states  of  the  same 
symmetry  at  the  same  energy  associated  with  the  majority-spin  electrons  upon  crossing  the 
Fermi  level.  We  obtain  a  band  offset  of  4.1  eV  for  the  A2  band  in  the  [100]  direction,  and 
another  band  offset  of  3.3  eV  for  the  Si  band  in  the  [110]  direction  [see  Fig.  6].  Ortega  et 
a/ [37]  have  shown  for  Co/Cu  and  Fe/Cu  systems  that  a  large  band  offset  between  the  band 
edges  of  the  magnetic  and  the  spacer  layer  gives  rise  to  sharply  confined  wavefunctions. 

For  2;  being  the  multilayer  growth  direction,  the  translational  symmetry  holds  only  in  the 
x-y  plane  so  that  the  wavefunction  of  an  electron  in  a  quantum  well  can  be  expressed  as [38] 

=  AMz)E^-Pu^,,,{R),  (6) 

where  A  is  a  normalization  constant;  k  =  {kx,ky)]  p  =  (x,r/);  Un^h{R)  is  periodic  in  the 
lattice  while  fk{z)  is  the  envelope  function,  which  ensures  that  the  boundary  conditions 
are  met  at  the  interfaces.  Eq.  (6)  shows  that  the  wavefunction  of  a  QW  state  consists  of  a 
rapidly  oscillating  Bloch  function,  which  is  modulated  by  a  slowly  varying  envelope  function. 
For  a  single  band  edge,  the  modulation  of  the  Bloch  wave  {k^dge)  by  an  envelope  {kenv)  yields 
a  total  wavevector[37] 

ktot  —  kgdge  i  ^envi  ('i^) 

where  only  the  normal  components  of  the  wavevectors  are  considered  since  the  boundary 
conditions  for  the  components  that  are  parallel  to  the  interface  mimic  the  bulk.  By  drawing 
an  analogy  with  the  simple  interferometer  model,  we  may  expect  the  interference  maxima 
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to  appear  at  every  half  wavelength  (Aenv)  of  the  envelope  as  the  thickness  of  the  Nb  spacer 
layer  is  increased  so  that  QW  states  appear  with  a  period  of  Aer^^;/2.  For  QW  states  at  the 
Fermi  level,  the  total  wavevector  becomes  the  Fermi  wavevector  if  the  spacer  layer  is  thick 
enough  to  exhibit  bulk-like  bonding[39].  As  Fig.  6(a)  shows,  the  band  edge  for  the  Si  band 
in  the  [110]  direction  is  located  at  the  zone  boundary  {ksz)  so  that  kgnv  =  {ksz  -  kp).  The 
frequency  associated  with  this  k^nv  in  real  space  may  be  viewed  as  a  beat  frequency  between 
the  Fermi  wavelength  and  the  atom  spacing[10].  The  oscillation  period  thus  becomes 

T  =  -K/kenv  =  7r/ {ksz  “  kp)  (8) 

This  way,  we  obtain  a  period  of  5.8  A  due  to  the  Si  band,  which  is  in  fairly  good  agreement 
with  the  periodicity  of  6.1  A,  as  appeared  in  the  oscillating  DOS  at  Ep  (see  Fig.  7).  It  may 
be  noted  that  the  wavevector  determining  the  periodicity  of  the  envelope  function  along  the 
[110]  direction  is  quite  small  and  thus  the  wavelength  is  large.  However,  a  reverse  situation 
occurs  for  the  A2  band  in  the  [100]  direction,  where  {kpz  -  kp)  is  large  as  kp  is  small.  The 
envelope  function  should  produce  in  this  case  a  short-wavelength  weak  modulation  of  the 
QW  wavevector  [38,  40]. 

E.  Phase  accumulation  model 

A  simple  way  of  predicting  the  energy  for  QW  states  is  the  phase  accumulation  model 
(PAM),  which  is  often  used  to  describe  the  thickness  dependence  of  QW  energies  and  reso¬ 
nances  in  thin  films  and  other  layered  structures [40].  According  to  this  model,  the  condition 
for  a  QW  state  is  determined  by  summing  over  all  the  phases  accumulated  by  a  propagat¬ 
ing  plane  wave  inside  a  quantum  well.  The  total  phase  accumulated  must  be  an  integral 
multiple  of  27r.  For  m  monolayers  of  Nb  spacer,  each  of  thickness  d,  the  total  spacer  width 
becomes  =  md.  Thus  the  distance  traversed  by  an  electron  in  the  spacer  medium  is 
2md  resulting  in  2mdk^  phase  change,  where  A;-'-  represents  the  electron  wavevector  normal 
to  the  layers.  If  the  phase  shift  of  the  electron  wavefunction  upon  reflection  at  each  interface 
is  (j)!,  the  total  phase  in  a  round  trip  within  the  spacer  medium  can  be  written  as[40] 

2k'^dpfb  +  ‘2(f)i{E)  =  2n7r,  (9) 

where  n  is  an  integer  related  to  the  number  of  half-wavelengths  that  span  the  quantum  well. 
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The  phase  shift  at  each  interface  in  Eq.  (9)  can  be  calculated  by  making  use  of  following 
ansatz[40] 

=  2  sin~^  ^ ~ 

where  Eu  and  E^  represent  the  energies  of  the  potential  well  taken  from  the  upper  and 
lower  edges  of  the  energy  gap  in  the  Fe  majority-spin  band,  which  acts  as  a  potential  energy 
barrier  for  the  propagating  electrons.  The  interface  reflection  phase  evolves  across  a  band 
gap  from  0  to  -tt  in  traversing  from  the  top  to  the  bottom  of  the  potential  well. 

In  order  to  obtain  A;-*-  in  Eq.  (9),  we  make  use  of  a  two-band  nearly  free  electron  (NEE) 
model,  which  approximates  the  Nb  electrons  in  the  A2  and  Si  bands.  According  to  this 
model,  the  wavevector  of  a  NEE  band  can  be  expressed  as[9,  41] 


where  G  =  where  m*  is  the  electron  effective  mass;  Vq  is  a  constant  offset  of  the 

periodic  potential;  2U  is  the  energy  gap  at  the  zone  boundary;  and  E  is  the  electron  energy 
with  respect  to  the  Fermi  level.  Eq.  (10)  thus  needs  three  parameters  to  determine  for  a 
given  energy  E.  A  fit  to  the  self-consistently  calculated  A2  band  of  Nb  yields  U  =  2.05  eV, 
Vo  =  -9.85  eV  and  m*  =  l.OSmg,  where  mg  is  the  electron  mass.  On  the  other  hand,  upon 
fitting  Eq.  (11)  to  the  Si  band  of  Nb,  we  have  m*  =  l.Obmg,  U  =  1.66  eV  and  Vq  =  -5.5 
eV. 

The  presence  of  a  periodic  atomic  potential  in  the  spacer  layer,  however,  ensures  that  in 
addition  to  two  traveling  waves  with  wavevectors  k  and  -k  arising  from  reflections  at  the 
interfaces,  electron  wavefunctions  in  the  phase  accumulation  model  should  be  described[40, 
42]  by  the  combination  of  two  more  waves  with  wavevectors  (k-g)  and  -{k-g),  where  g  is  the 
reciprocal  lattice  vector.  These  additional  waves  correspond  to  the  Bragg  scattering  within 
the  periodic  potential  of  the  spacer  layer.  Since  the  reciprocal  lattice  vector  normal  to  the 
interface  within  the  hrst  Brillouin  zone  is  g  =  2kBz  =  2'K/d,  we  have  2kBzdNb  =  2m7r.  Using 
this  relation,  Eq.  (9)  can  be  rewritten  as 

2K“‘-div6  —  2(f)i{E)  =  2i'7t,  (12) 

where  =  kBz  -  and  u  =  m  -  n.  It  follows  from  Eq.  (8)  of  the  preceding  section  that 
K-^  is  characteristic  of  the  envelope  function  that  modulates  the  QW  wavefunction. 
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As  we  know,  states  within  Nb  layers  that  do  not  coincide  in  energy  and  momentum 
with  those  within  Fe  layers  form  the  quantum  well  states.  With  the  increase  in  the  Nb 
layer  thickness,  the  positions  of  the  spin  polarized  QW  states  vary  and  thus  exhibit  regular 
dispersion  through  the  Fermi  level.  For  n  =  1,  Figs.  8(a)  and  8(b)  demonstrate  the  nature 
of  QW  dispersions  in  Fe/Nb  multilayers  along  [100]  and  [110]  directions  as  generated  by  Eq. 
(12)  of  the  phase  accumulation  model.  We  find  that  the  A2  band  yields  weak  dispersion 
as  compared  to  the  Si  band  of  Nb  due  mainly  to  the  weak  modulation  of  A;-*-  in  the  [100] 
direction.  The  oscillation  periods  of  the  quantum  well  states  at  the  Fermi  energy,  as  shown 
by  the  dispersion  curves  of  8(a)  and  8(b),  are  obtained  as  about  5.5  and  7.3  A.  These  values 
are  close  to  the  fitted  values  of  6.3  and  7.7  A  respectively. 

IV.  CONCLUSIONS 

In  this  paper,  we  have  explained  the  phenomena  of  interlayer  exchange  coupling  in 
Fe/Nb (001)  multilayers  in  terms  of  RKKY  as  well  as  the  QW  model  using  the  density 
functional  calculations.  The  RKKY  periodicities  arising  out  of  extremal  spanning  vectors 
of  the  multisheet  Nb  Fermi  surface  have  been  found  to  be  in  favorable  agreement  with  the 
available  experimental  as  well  as  calculated  data.  In  the  region  of  small  magnetic  moments  of 
Fe,  both  the  short  and  long  periods  display  the  RKKY  kind  of  coupling,  since  lEC  shows  bi¬ 
linearity  in  magnetization.  However,  the  presence  of  additional  biquadratic  and  triquadratic 
coupling  constants  in  lEC  at  higher  Fe  moments  signify  the  onset  of  the  non-RKKY  char¬ 
acter,  especially  of  the  short-period  oscillations.  We  have  found  that  at  moments  closer 
to  the  bulk  value,  the  long-period  oscillations  become  eventually  independent  of  the  size 
of  the  magnetic  moments.  The  appearance  of  higher  harmonics  in  the  well-fitted  envelope 
of  the  oscillatory  exchange  coupling  turns  out  to  be  responsible  for  the  saturation  of  long 
periods.  The  oscillatory  nature  of  the  density  of  states  at  the  Fermi  level  supports  the  QW 
description  of  the  exchange  coupling  in  Fe/Nb  multilayers.  It  is  because  such  oscillations 
originate  mainly  from  the  spin-dependent  confinement  of  the  propagating  electrons  inside 
the  quantum  well  of  the  spacer  medium.  Quantum  well  dispersions  around  the  Fermi  level 
illustrate  that  the  majority  spin  bands  contribute  largely  to  the  formation  of  QW  states, 
which  is  subsequently  analyzed  by  making  use  of  the  phase  accumulation  model.  All  these 
results  show  that  magnetic  quantum  wells  are  formed  in  Fe/Nb  multilayers  and  the  oscil- 
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latory  behavior  of  the  exchange  coupling  in  Fe/Nb  multilayers  is  better  described  by  the 
Quantum  Well  model. 
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Figure  Captions 


FIG.  1:  The  calculated  oscillatory  interlayer  exchange  coupling  (solid  circles)  as  a  function  of  the 
number  of  Nb  spacer  layers  in  Fe3NbTO(m  =  1-16)  multilayers.  The  solid  line  is  the  fitted  plot  (see 
the  text  for  details). The  thickness  of  1  Nb  layer=1.5335  A. 

FIG.  2:  (Color  online)  Cross  sections  of  the  Fermi  surface  of  Nb  in  the  (100)  plane.  F  labels  the 
center  of  the  Brillouin  zone,  N  indicates  the  center  of  each  face  of  the  dodecahedron  and  H  labels 
the  corners  of  the  four-fold  symmetry  on  the  zone  boundary. 

FIG.  3:  Variation  of  the  Fe  magnetic  moment  with  Nb  spacer  thickness  in  Fe3Nb^(m  =  1-16) 
heterostructures  (see  the  text  for  details).  The  inset  shows  the  experimental  results  by  Mattson  et 
al[13]. 

FIG.  4:  The  induced  magnetic  moment  in  Nb  spacer  layers  for  Fe3Nbi6  heterostructures  in  the 
(a)  ferromagnetic  (FM)  and  (b)  antiferromagnetic  (AFM)  orientations.  The  thickness  of  1  Nb 
layer=1.5335  A. 

FIG.  5:  (Color  online)  Interlayer  exchange  coupling  in  Fe/Nb  multilayers  as  a  function  of  Fe 
magnetic  moment  parametrized  by  a  (see  the  text  for  details). The  asterisks  represent  the  results 
for  the  16-atom  supercell  while  the  solid  line  depicts  the  fitted  curve  according  to  Eq.  (5).  Note 
that  the  linear  part  of  the  solid  line  resembles  the  RKKY  kind  of  coupling.  The  triangles  and 
diamonds  indicate  the  short-period  amplitudes  while  the  circles  and  squares  represent  the  long- 
period  amplitudes  of  the  Fe3Nb^(m  =  1-16)  multilayer  system. 

FIG.  6:  The  bulk  energy  bands  of  Nb  and  bcc  Fe  (t  and  4-)  along  the  (a)  [100]  and  (b)  [110] 
directions.  The  bands  with  A2  and  Si  symmetries  are  displayed  by  solid  lines  while  the  bands 
with  Si  symmetry  by  dashed  lines.  Only  majority-spin  states  in  Nb  exhibit  quantum  well  character 
at  the  Fermi  level  since  the  minority-spin  A2  and  Si  states  couple  with  the  corresponding  states 
in  Fe. 

FIG.  7:  Oscillations  in  the  density  of  states  at  the  Fermi  level,  Ep^  with  the  Nb  spacer  thickness, 
caused  by  the  quantum  well  states  in  Fe3NbTO(m  =  1-16)  heterostructures. 

FIG.  8:  The  thickness  dependence  of  the  QW  energies  in  Fe/Nb  multilayers  generated  by  Eq.  (12) 
of  the  phase  accumulation  model  with  respect  to  (a)  A2  and  (b)  Si  bands. 
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